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ABSTRACT 

An inhomogeneous cosmological magnetic field creates vortical perturbations that 
survive Silk damping on much smaller scales than compressional modes. This ensures 
that there is no sharp cut-off in anisotropy on arc-minute scales. As we had pointed 
out earlier, tangled magnetic fields, if they exist, will then be a potentially impor- 
tant contributor to small-angular scale CMBR anisotropies. Several ongoing and new 
experiments, are expected to probe the very small angular scales, corresponding to 
multipoles with I > 1000. In view of this observational focus, we revisit the predicted 
signals due to primordial tangled magnetic fields, for different spectra and different 
cosmological parameters. We also identify a new regime, where the photon mean-free 
path exceeds the scale of the perturbation, which dominates the predicted signal at 
very high I. A scale- invariant spectrum of tangled fields which redshifts to a present 
value B = 3 x 10~ 9 Gauss, produces temperature anisotropies at the 10/iK level 
between I ~ 1000 — 3000. Larger signals result if the univese is lambda dominated, 
if the baryon density is larger, or if the spectral index of magnetic tangles is steeper, 
n > —3. The signal will also have non-Gaussian statistics. We predict the distinctive 
form of the increased power expected in the microwave background at high / in the 
presence of significant tangled magnetic fields. We may be on the verge of detecting or 
ruling out the presence of tangled magnetic fields which are strong enough to influence 
the formation of large-scale structure in the Universe. 
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1 INTRODUCTION 

The origin of large-scale cosmic magnetic fields remains a 
challenging problem. It is widely assumed that magnetic 
fields in astronomical objects, like galaxies, grew by tur- 
bulent dynamo action on small seed magnetic fields (cf. 
Ruzmaikin, Shukurov & Sokoloff 1988; Beck et al 1996). 
However, the efficiency of turbulent galactic dynamos is still 
unclear, especially in view of the constraints implied by he- 
licity conservation (Cattaneo & Vainshtein 1991; Kulsrud 
& Anderson 1992; Gruzinov & Diamond 1994; Cattaneo & 
Hughes 1996;Subramanian 1998, 1999; Blackman & Field 
2000; Kleoorin et al 2000; Brandenburg 2001; Brandenburg 
and Subramanian 2000; Brandenburg, Dobler and Subrama- 
nian 2002). Magnetic fields with larger coherence scales may 
also be present in clusters of galaxies (Clarke, Kronberg & 
Bohringer 2001) and at high redshifts (Oren & Wolfe 1995). 
It is not yet clear that all these large-scale coherent fields 
could result from dynamo action. Alternatively, galactic or 
cluster fields could be remnants of a primordial cosmologi- 



cal magnetic field (cf. Kulsrud 1990; 1999), although, as yet, 
there is no entirely compelling mechanism for producing the 
required field. It could be present in the initial conditions, 
be produced quantum gravitationally or at a phase transi- 
tion, or be generated in some way at the end of a period of 
inflation, perhaps with an almost scale-invariant spectrum 
(Turner & Widrow 1988; Ratra 1992; cf. Grasso and Ruben- 
stein 2001 for a review). 

A primordial field that expanded to contribute a present 
field strength of order 10 -9 Gauss, tangled on galactic scales, 
could also affect the process of galaxy formation (Rees & 
Reinhardt 1972; Wasserman 1978; Kim, Olinto & Rosner 
1996; Subramanian & Barrow 1998a, SB98 hereafter). It is 
of considerable interest, therefore, to find different ways of 
limiting or detecting such primordial fields (see Kronberg 
1994 and Grasso & Rubenstein 2001 for reviews). Here, we 
will show that the imminent extension of cosmic microwave 
background radiation (CMBR) observations to very small 
angular scales allows us to test for the presence of dynam- 
ically significant cosmological magnetic fields in new ways. 
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This is because tangled magnetic fields create small-scale 
power in CMBR temperature anisotropy that would have 
been damped out if magnetic fields were absent. Conversely, 
if enhanced power is found in the CMBR power spectrum on 
very small scales then it might be explained by the effects 
of tangled magnetic fields. 

In an earlier paper (Subramanian & Barrow 1998b, 
hereafter Paper I), we argued that observations of 
anisotropies in the CMBR provide a potentially power- 
ful constraint on such tangled magnetic fields. Indeed, the 
isotropy of the CMBR already places a limit of 6.8 x 
10~ 9 (Q m h 2 ) 5 Gauss on the present strength of any uniform 
(spatially homogeneous) component of the magnetic field 
(Barrow, Ferreira & Silk 1997), where Q m is the present 
matter density parameter, and h the Hubble constant in 
units of 100 km s _1 Mpc -1 ). In Paper I we obtained com- 
parable constraints on tangled (inhomogeneous) magnetic 
fields and highlighted the distinctive fluctuation signature 
that they are expected to leave in the small-scale structure 
of the CMBR. In particular, tangled magnetic fields pro- 
duce vortical perturbations, which are overdamped in the 
radiation era. These can then survive Silk damping (Silk 
1968) on scales much smaller than the compressional modes 
(Jedamzik, Katalinic & Olinto 1998; SB98), so an obvious 
place to expect signals induced by tangled magnetic fields 
is below the Silk damping scale, or at multipoles of greater 
than I ~ 1000. Several new and ongoing experiments (MAP, 
VSA, ACBAR, CBI, ATCA and Planck Surveyor) are indeed 
expected to provide information in this large I regime. This 
motivates us to revisit the computation of the expected sig- 
nals at larger values of I, for a wider variety of cosmological 
parameters and spectral indices, than were made in Paper 
I. In particular, we calculate results for a new high-Z regime, 
where the photon mean-free path exceeds the scale of the 
perturbation. These predictions will allow comparison with 
future observations. 



2 ESTIMATES OF THE INDUCED 
ANISOTROPY 

The evolution of temperature anisotropy for vector pertur- 
bations has been derived in detail by Hu and White (1997a) 
in the total angular momentum representation and was also 
given in Paper I. The temperature anisotropy is expanded in 
terms of tensor spherical harmonics and the angular power 
spectrum induced by vector perturbations is given by (see 
Paper I and Mack, Kahniashvilli and Kosowsky 2002), 



Ci = 4tt 



x < 



k 2 dk 1(1 + 1) 
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Here, v(k, r) is the magnitude of the vorticity, Qi 



v, 



in Fourier space, where vf is the rotational component of 
the fluid velocity, and Vi is the vector metric perturbation. 
Note that since ili also appears in the Euler equation for 
baryons, we need never compute the vector metric perturba- 
tion explicitly (see also Mack, Kahniashvilli and Kosowsky 
2002). Also, k is the co-moving wave number, r is conformal 
time, to its present value, and ji(z) is the spherical Bessel 
function of order I. We have ignored a small polarization 



correction to the source term and also a metric perturba- 
tion term which is subdominant at large I (cf. Seshadri and 
Subramanian 2001). The 'visibility function', g(ro, r), deter- 
mines the probability that a photon reaches us at epoch to if 
it was last scattered at the epoch r. We adopt a flat universe 
throughout, with a total matter density tt m and a non-zero 
cosmological constant density Q,a = 1 — fim today. 

In Paper I we obtained analytic estimates for C'i in var- 
ious asymptotic regimes. We briefly recapitulate the argu- 
ments and results. Firstly, we approximated the visibility 
function as a Gaussian: <?(to,t) = (2-ko 2 )~ x ^ 2 exp[— (t — 
T») 2 /(2cr 2 )], where r* is the conformal epoch of "last scatter- 
ing" and a measures the width of the LSS. To estimate these, 
we use the expressions given in Hu and Sugiyama (1995), 
adopting a baryon density parameter fif, = 0.02/i~ 2 , and h — 
0.7. We estimate the epoch of last scattering, to be (l + z») ~ 
1150 and a/r» ~ 0.064. To convert redshift into conformal 
time we use r — 6000/t -1 ((a + a e9 ) 1/ ' 2 — a]/ q 2 ) /qU. 2 , valid for 
a flat universe (cf. Hu & White 1997b). Here, the expansion 
factor a = (l + z) -1 and a eq = 4.17 x 10~ 5 (fi m /i 2 )" 1 . For an 
fi m = 1 model, we get r» = 130.0/ir 1 Mpc, and a = 8.4/tT 1 
Mpc, while for a A-dominated model with £l m = 0.3, we get 
t* = 187.5/i _1 Mpc and a = 12.1/i _1 Mpc. We will use these 
numbers in the numerical estimates below. 

The dominant contributions to the integral over t in 
Eq. (|l|) then come from a range a around the epoch r = t* . 
Furthermore, ji(k(ro — t)) picks out (k,r) values in the in- 
tegrand which have k(ro — t) ~ I. Thus, following the ar- 
guments detailed in Paper I, for ka « 1 we get the an- 
alytical estimate, 1(1 + 1)C;/(2tt) « (tt/4) A 2 (k, T*)\ k =i/R. ■ 
Here, A 2 — k s < |t>(fc, t»)| 2 > /(2ty 2 ) is the power per unit 
logarithmic interval of k, residing in the net vorticity per- 
turbation, and R* — to — r*. In the opposite limit, ka >> 1, 
we get 1(1 + l)C,/(2n) « (0F/4)(A 2 (fc, T*)/(ka)\ k=l/R , . At 
small wavelengths, C; is suppressed by a 1/ka factor due to 
the finite thickness of the LSS. 

To evaluate Ci, one also needs to estimate v, the vor- 
ticity induced by magnetic inhomogeneities. We assume the 
magnetic field to be initially a Gaussian random field. On 
galactic scales and above, the induced velocity is gener- 
ally so small that it does not lead to any appreciable dis- 
tortion of the initial field (Jedamzik, Katalinic and Olinto 
1998, SB98). So, to a very good approximation, the mag- 
netic field simply redshifts away as B(x,t) = bo(x)/a 2 . 
The Lorentz force associated with the tangled field is then 
Fl = (V x bo) x bo/(4-7ra 5 ), which pushes the fluid and 
creates rotational velocity perturbations. These can be es- 
timated as in Paper I, by using the Euler equation for the 
baryons. On scales larger than the photon mean-free-path 
at decoupling, where the viscous effect due to photons can 
be treated in the diffusion approximation, this reads 



4 \ dfli 



p\y da 
a dt 



k 2 n 



p. . p. 

Q _ 3 

47ra b 



(2) 



Here, p 7 is the photon density, pb the baryon density, 
and rj — (4/15)p 7 Z 7 the shear viscosity coefficient asso- 
ciated with the damping due to photons, whose mean- 
free-path is Z 7 = (rigor) 1 = L y a(t), where n e is the 
electron density and ot the Thomson cross-section. For 
(1 + z*) ~ 1150, we get L 7 (t») ~ 1.8/ 6 _1 Mpc, where 
fb = (0.bh 2 /0.02). We have defined the Fourier transforms 
of the magnetic field, by bo(x) = ]T] k b(k) exp(ik.x) and 
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F(k) = E p [b(k + p).b*(p)]p-[k.b*(p)]b(k + p). The pro- 
jection tensor, Pij(k) = [5ij — kikj/k 2 ] projects F onto its 
transverse components perpendicular to k. 

The comoving Silk damping scale at recombination, 
Ls = kg 1 ~ 10 Mpc, separates scales on which the ra- 
diative viscosity is important (kLs 3> 1) from those on 
which it is negligible (kLs <§C 1)- For kL s <C 1, the damp- 
ing due to the photon viscosity can be neglected compared 
to the Lorentz force. Integrating the baryon Euler equa- 
tion, assuming negligible initial vorticity perturbation, then 
gives Qi = Gir/(1 + S*), where d = 3PijFj /[Wirpo], po 
is the present-day value of p 7 , and S* = (3pft/4p 7 )(r*) ~ 
0.53(f2 6 ft 2 /0.02). In the other limit, with kL 3 >> 1, we can 
use the terminal-velocity approximation, neglecting the in- 
ertial terms in the Euler equation, to balance the Lorentz 
force by friction. This gives Qi = (Gi/k)(5/kL 7 ), on scales 
where diffusion damping operates. The transition Silk scale 
can also be estimated by equating Qi in the two cases, to 
give fc s ~ [5(1 + S*)/(tL 1 (t))] 1 ' 2 . 

A new regime arises on very small scales that are well 
below the photon mean-free path. The radiative drag force 
is then no longer described by the diffusion approximation, 
but rather by the free-streaming drag given in Eq. (6.1) 
of SB98. Under strong damping, this term dominates the 
inertial terms in the Euler equation and the fluid reaches 
terminal velocity where friction due to free-streaming pho- 
tons balances driving by the Lorentz force, Fl- From Eq. 
(6.2) of SB98, we get 4p 7 / (3p b )n e a T £l = F L /p b , which 
gives Qi = GiL-y. One can also estimate the transition scale, 
say kjl, below which free-streaming damping dominates, 
by equating the induced velocities in the diffusion and free- 
streaming damping regimes, to get kf s ~ VE/L-y. 

In order to compute the Cis we also need to specify 
the spectrum of the tangled magnetic field, say M(k). We 
define, < 6i(k)6j(q) >= 5k,q-Py (k)M(fe), where <5k, q is the 
Kronecker delta which is non-zero only for k = q. This gives 
< b 2 , >= 2f(dk/k)Al(k), where A 2 (fe) = k 3 M(k)/(2n 2 ) 
is the power per logarithmic interval in k space residing 
in magnetic tangles, and we replace the summation over 
k space by an integration. The ensemble average < \v\ 2 >, 
and hence the C;s, can be computed in terms of the mag- 
netic spectrum M(k). It is convenient to define a dimension- 
less spectrum, h(k) = A 2 .(k)/(B 2 /2), where Bo is a fiducial 
constant magnetic field. The Alfven velocity, Va, for this 
fiducial field is, 



Bo 



3.8 x 10" 4 £- 



(3) 



(167TP0/3) 1 / 2 

where B-g = (_Bo/lCP 9 Gauss). We will also consider power- 
law magnetic spectra, M(k) = Ak n cut-off at k — k c , where 
fc c is the Alfven-wave damping length-scale (Jedamzik, 
Katalinic and Olinto, SB98). We fix A by demanding that 
the smoothed field strength over a "galactic" scale, ka = 
lftMpc -1 , (using a sharp fe-space filter) is Bo, giving a di- 
mensionless spectrum for n > — 3 of 



h(k) = (n + 3)(k/k G ) 3+n . 



(4) 



We can now put together the above results to derive an- 
alytic estimates for the CMBR anisotropy induced by tan- 
gled magnetic fields. As a measure of the anisotropy we 
define the quantity AT (I) = [1(1 + 1)Ci/2tt] 1/2 T , where 
T = 2.728 K is the CMBR temperature. On large scales, 



such that kL s < 1 and ka < 1, the resulting CMBR 
anisotropy is (see Paper I) 

AT B (i) = To(g)V»I(fc )7 *^ 



5.8pK 



B- 



(1 + 5,) 



I 



(•») 



Here, I = kR* and we have used cosmological parameters 
for the A-dominated model, with Qa = 0.7, Q. m = 0.3 and 
Q. b h 2 = 0.02 (in Paper I, we used a purely matter-dominated 
fim = 1 model). We also use the fit given by Hu and White 
(1997b) to calculate r = 6000ft- 1 ((1 + a eq ) 1/2 - a^ 2 )(l - 
0.0841Zn(!T2 m ))/^ /2 , valid for flat universe. 

On scales where kLs > 1 and ka > 1, but fcL 7 (r») < 1, 
we get 

7T 1/4 W 2 

AT B (l) = T ?-=I(k) ur ,\f. , 1/a 
^32 y fcL 7 (r,)(fccr) 1 /2 



B. 



I 



2000 



-3/2 



where h-;o = (ft/0.7) and ft = (Qth /0.02); there is also 
a weaker dependence of other parameters on /(, but the 

oc ft 1 - 



strongest dependence comes from the fact that L 7 
In the free-streaming regime, the estimate for AT (I) is ob- 
tained by replacing (feL 7 /5) _1 in Eq. (^|) by (fcL 7 ). The 
CMBR anisotropy which results for scales so small that free- 
streaming damping dominates is 



ATfl(i) = 



V32 (ka) 1 / 2 



(r.) 



0.4^ 



B-g 
3 



/ 



20000 



1/2 



-T- (7) 



The function I 2 (k) in the Eqs.< 
mode-coupling integral given by 



-((?]) is a dimensionless 



I 2 (k) 



dq 

q 

x(l-M 2 ) 



dp, 



ft(g)ft(|(k + q)|)fc 3 



(fe 2 - 

(k 



q 2 + 2kqp) 3 / 2 
- 2qp)(k + qp) 



(k 2 + q 2 + 2kq[i) 



where |(k+q| = (k 2 + q 2 +2kqp) 1 ^ 2 . In general, I(k) can only 
be evaluated numerically but for h(k) = k5o(k — ko), where 
5d(x) is the Dirac delta function, it can be evaluated exactly. 
One gets, I(k) = (k/k )[l - (k/2k ) 2 ] 1/2 , for k < 2k , and 
zero for larger k. So in this case, I(k) contributes a factor of 
order unity around k ~ ko, with I(ko) = v3/2. 

We can also find an analytic approximation to the mode 
coupling integral for power-law magnetic spectra. The ap- 
proximation is different for n > —3/2, and for n < —3/2. For 
n > —3/2 and for k « k c (which is relevant for / << k c R,), 
one gets (Seshadri & Subramanian 2001): 



I 2 (k) 



28 (n + 3) 2 . k 3 k c 3+2n 
15 (3 + 2n) ( k G k G 



(9) 



The mode-coupling integral is dominated by the small scale 
cut-off in this case. In the other limit n < —3/2, we find: 

J 2 (fc) = |(n + 3)(A)6+ 2 " (io) 

For a nearly scale invariant spectrum, say with n = —2.9, 
we then get AT (I) ~ 4.7/iJi'(;/f000) 1 ' 1 for scales larger than 
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the Silk scale, and AT {I) ~ 5.6m#(//2000) _1 ' 4 , for scales 
smaller than Ls but larger than L 7 . Larger signals will be 
expected for steeper spectra, n > —2.9 at the higher / end. 

To complement these analytic results, we have also com- 
puted AT (I) for the above spectra, by evaluating the r 
and k integrals in Eq.([jJ) numerically. We retain the ana- 
lytic approximations to I(k) and fij(fc), with transitions be- 
tween the limiting forms at wavenumbers ks and kf 3 . The 
results are displayed in Figures 1 and 2. We see that for 
Bo ~ 3 x 10 -9 G, this leads to a predicted RMS temperature 
anisotropy in the CMBR of order 10/j.K for 1000 < I < 3000, 
for a nearly scale-invariant power law spectra with n = —2.9. 
Larger signals result, at the high I end, for a A-dominated 
universe, compared to a matter dominated model (compare 
the solid and dotted curves), basically due to an increase in 
i?* for this model. Larger signals also result for larger baryon 
density (compare the solid and dashed-triple-dotted curves), 
due to the decrease in L 7 and hence the damping effects of 
radiative viscosity. Also, a moderately steeper spectral in- 
dex, with n — —2.5, which has more power on small scales 
leads to an increased AT (compare the dashed-triple-dotted 
and dashed curves). Much larger signals result from even 
steeper spectra, but we have not displayed these results, as 
spectra with n > —2.5 and B-g ~ 3, are probably ruled 
out because of gravitational wave production, estimated by 
Caprini and Durrer (2002). They will also be significantly 
constrained by the high-/ limits on the anisotropy by ATCA, 
given in Subrahmanyan et al (1998). We find that our ana- 
lytic approximation of the r and k integrals in Eq.(|l|), tends 
to underestimate the amplitude of AT (I) by about a fac- 
tor ~ 2, although the analytically predicted /-dependences 
agree very well with the numerical integration at both small 
and large I. The numerical integration over the Bessel func- 
tion and the visibility functions are of course expected to be 
more accurate. The new regime of free-streaming damping 
is only seen at very large I ~ kf s R„, which is of order 10 4 
even for the matter-dominated model, and is ~ 20, 000 for 
the A-dominated model. 



3 DISCUSSION 

We re-examined the expected small-angular scale anisotropy 
induced by tangled magnetic fields, for different cosmolog- 
ical parameters and spectral indices and including a new 
regime at very small scales. We are motivated primarily by 
several ongoing and future experiments, which probe this 
large- I regime. Tangled magnetic fields give a distinctive 
contribution to the high-/ signal, since they create vortical 
perturbations that survive Silk damping on much smaller 
scales than do compressional modes. Moreover, for small- 
scale rotational perturbations, the damping due to the fi- 
nite thickness of the LSS is also milder than for compres- 
sional modes. It should be emphasized that, by contrast, 
in the standard non- magnetic models the C;s have a sharp 
cut off for / > ksR*, due to Silk damping. As we see from 
Figures 1 and 2, it is precisely here that the magnetically 
induced signals begin to dominate. A scale-invariant spec- 
trum of tangled fields which redshifts to a present value 
Bo = 3 X 10 -9 Gauss, produces temperature anisotropies 
at the 10/j.K level between I ~ 1000 — 3000. Larger signals 
are produced in a A-dominated universe, in a universe with 
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Figure 1. AT versus / predictions for different cosmological mod- 
els and magnetic power spectrum M(k) oc k n , for B_g = 3. The 
bold solid line ( — ) is for a canonical fiat, A-dominated model, 
with f2 A = 0.7, n m = 0.3, tt b h 2 = 0.02, h = 0.7 and almost 
scale invariant spectrum n = —2.9. The dotted curve (....) ob- 
tains when one changes to Q m = 1 and Qa = model. The 
dashed-triple-dotted curve (-■ ■ — ) is for the A-dominated model 
with a larger baryon density Q^h 2 = 0.03, while the dashed curve 

( ) changes this model by adopting a magnetic spectral index 

of n = —2.5. We also show for qualitative comparison (dashed- 
dotted curve — • — ), the temperature anisotropy in a 'standard' 
A-CDM model, computed using CMBFAST (Seljak & Zaldar- 
riaga 1996), with cosmological parameters as for the first model 
described above. These curves show the build up of power due to 
vortical perturbations from tangled magnetic fields which survive 
Silk damping at high I ~ 1000 — 3000. The eventual slow decline 
is due to the damping by photon viscosity, although this decline 
is only a mild decline as the magnetically sourced vortical mode 
is overdamped. By contrast, in the absence of magnetic tangles 
there is a sharp cut off due to Silk damping. 



larger baryon density, or if the spectral index of magnetic 
tangles is steeper, n > —3. Note that the anisotropy in hot 
or cold spots could be several times larger, because the non- 
linear dependence of Ci on M(k) will imply non-Gaussian 
statistics for the anisotropies. Indeed, as we pointed out in 
Paper I , the AT (I) will obey a modified form of the % 2 
distribution, the probability falling off almost exponentially, 
rather than as a Gaussian. 

We see that the predicted signal is quite sensitive to 
both the cosmology and also the baryonic content (the value 
and any clumping of Q.b ) • It is also necessary in future work 
to calculate the damping effects more accurately, and also 
to take account of the effects of non-Gaussianity. The latter 
point is especially important in searches which cover lim- 
ited regions of the sky, and the variance may not be the 
best statistic to compare with observations. It is likely that 
the interpretation of the data at these small angular scales 
will also be complicated by the need to understand the con- 
tribution from discrete sources. 
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Figure 2. Log(AT) versus Log(l) predictions up to very large 
I ~ 10 4 ; the line codings are same as in Figure 1. 

The other potentially important source of temperature 
anisotropics at these small-angular scales is the Sunyaev- 
Zeldovich signal. However, it should be possible to isolate 
this signal by its frequency dependence. An important means 
of distinguishing the effects of tangled magnetic fields is by 
looking at the polarisation of the CMBR. it was shown by 
Seshadri and Subramanian (200f) that the vorticity induced 
by tangled magnetic fields produces significant and predom- 
inantly B-type polarization, at the micro-Kelvin levels, and 
this will also be dominant at large I. These should therefore 
be distinguishable from anisotropics produced by inflation- 
ary scalar and tensor perturbations. We hope to revisit this 
issue in a future study. 

Although we have concentrated on the vorticity induced 
anisotropy which is dominant on small-angular scales, tan- 
gled magnetic fields also produce anisotropies on larger an- 
gular scale, dominated by tensor metric perturbations in- 
duced by anisotropic magnetic stresses (Durrer, Ferreira & 
Kahniashvili 2000; Mack, Kahniashvili & Kosowsky 2002). 
Using the formalism described in these papers, we estimate 
a tensor contribution at small I < f00 of ~ 10.9/iK (Z/100) 01 
for the nearly scale invariant spectrum with n = —2.9 and 
AT (I) ~ 4.9fj,K(l/WQ) 0B for n = -2.5 and B_ 9 ~ 3. Since 
we have to add this power to the standard power produced 
by inflationary scalar perturbations, in quadrature, a tan- 
gled field with B_9 ~ 3 will produce of order a few to f0 
percent perturbation to the power in the standard CMBR 
anisotropy at small is. So if they are indeed detected at large 
I, below the Silk damping scale, and if they have a power- 
law spectrum extending to large-scales, one will also have to 
consider their effects seriously at large-angular scales, espe- 
cially in cosmological parameter estimation. 

We also note that a magnetic field which redshifts to a 
present day value of order a nano Gauss could impact sig- 
nificantly on the formation of galaxies (Rees & Reinhardt 
1972, Wasserman 1978, SB98). Our results therefore pro- 



vide a way of detecting the presence of small-scale magnetic 
inhomogeneities at a level which affects the formation of 
galaxies and clusters. 
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